Dispersion management for atomic matter waves 
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We demonstrate the control of the dispersion of matter wave packets utilizing periodic potentials. 
This is analogous to the technique of dispersion management known in photon optics. Matter wave 
packets are realized by Bose-Einstein condensates of ^'^Rb in an optical dipole potential acting as a 
one-dimensional waveguide. A weak optical lattice is used to control the dispersion relation of the 
matter waves during the propagation of the wave packets. The dynamics are observed in position 
space and interpreted using the concept of effective mass. By switching from positive to negative 
effective mass, the dynamics can be reversed. The breakdown of the approximation of constant, as 
well as experimental signatures of an infinite effective mass are studied. 



The broadening of particle wave packets due to the free 
space dispersion is one of the most prominent quantum 
phenomena cited in almost every textbook of quantum 
mechanics. The realization of Bose-Einstein condensates 
of dilute gases allows for the direct observation of wave 
packet dynamics in real space on a macroscopic scale Q . 
Using periodic potentials it becomes feasible to experi- 
mentally study to what extent the matter wave disper- 
sion relation can be engineered. This approach is similar 
to dispersion management for light pulses in spatially pe- 
riodic refractive index structures 2] . 

First experiments in this direction have already been 
undertaken in the context of Bloch oscillations of thermal 
atoms and condensates 0. The modification of the 
dipole mode oscillation frequency of a condensate due 
to the changed dispersion relation in the presence of a 
periodic potential has been studied in detail 0, 01- In 
contrast to these experiments where the center of mass 
motion was studied, we are investigating the evolution 
of the spatial distribution of the atomic cloud in a quasi 
one-dimensional situation. For the wave packet dynam- 
ics many new effects are expected due to the interplay 
between nonlinearity resulting from atom-atom interac- 
tion and dispersion. In particular, non-spreading wave 
ackets such as gap solitons and self trapped states 
are predicted. Our experiments show that even prop- 
agation in the linear regime with modified dispersion can 
change the shape of an initially Gaussian wave packet 
significantly. 

For atomic matter waves inside a one-dimensional opti- 
cal waveguide, we have achieved dispersion management 
by applying a weak periodic potential with adjustable ve- 
locity. Fig. n] shows the results of an experiment in which 
the propagation of an atomic wave packet is studied in 
the normal (Fig. ^p) and anomalous (Fig. ^) dispersion 
regime corresponding to positive and negative effective 
mass, respectively. A broadening of the wave packet is 
observed in both cases. The faster spreading in the case 
of anomalous dispersion is a consequence of the smaller 
absolute value of the negative effective mass. However, if 
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FIG. f : Controlling the dispersion of an atomic wave packet 
in a waveguide using a periodic potential. Shown are absorp- 
tion images of the wave packet averaged over four realizations 
(left) and the corresponding density distributions n(a;, t) along 
the waveguide (right), (a) Initial wave packet. (b,c) Images 
taken after an overall propagation time oi t = 26 ms for dif- 
ferent dispersion regimes with different effective masses as 
indicated, (d) Wave packet subjected to dispersion manage- 
ment: an initial stage of expansion for t = 17 ms with normal 
dispersion is followed by propagation with anomalous disper- 
sion for t = 9 ms. The broadening in the normal dispersion 
regime has been reversed by anomalous dispersion. 



we switch from one regime to the other during the propa- 
gation by changing the velocity of the periodic potential, 
the effects of normal and anomalous dispersion cancel. 
The wave packet, which has initially broadened under 
the influence of normal dispersion, reverses its expansion 
and compresses until it regains its initial shape (Fig.^). 
This is a direct proof of the realization of negative effec- 
tive mass. 

The concept of effective mass TOoff allows to describe 
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the dynamics of matter wave packets inside a periodic 
potential in a simple way via a modified dispersion rela- 
tion. The periodic potential in our experiments is well 
described by 



V{x) 



cos [Gx) 



with a modulation depth Vq on the order of the grating 
recoil energy Eq — fi^G^/d>m^ with G = 27r/(i where 
d — 417 nm represents the spatial period. The energy 
spectrum of atoms inside the periodic potential exhibits 
a band structure En{q) which is a periodic function of 
quasimomentum q with periodicity G corresponding to 
the width of the first Brillouin zone (Fig. In our 

experiment, we prepare condensates in the lowest energy 
band (n = 0) with a quasimomentum distribution 
centered aX q — with an rms width Aq <C G [loj . 

It has been shown by Steel et al. that in 

this case the condensate wave function in a quasi one- 
dimensional situation can be described by = 
A(x,t) (l)q^{x) exp{—iE{)[qc)t/K), where represents 
the Bloch function in the lowest energy band correspond- 
ing to the central quasimomentum. The evolution of the 
envelope function t), normalized to the total number 
of atoms A'o, is described by 
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The strength of the atom-atom interaction is given by 
g = o:„i2hLL}±a^ with the transverse trapping frequency 
of the waveguide uj±, the s-wave scattering length a, 
and a renormalization factor a„; = l/d dx\(f)q^\^ . Be- 
sides the modification of the nonlinear term the peri- 
odic potential leads to a group velocity of the enve- 
lope A{x, t) determined by the energy band via Vg{qc) — 
h~^{dEa{q)/dq)\q^ (Fig.Eb). In addition, the kinetic en- 
ergy term describing the dispersion of the wave packet is 
modified by the effective mass (Fig. ^) 
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It is important to note that group velocity and effective 
mass are assumed to be constant during the propagation. 
Since the approximation of constant effective mass corre- 
sponds to a parabolic approximation of the energy band, 
it is only valid for sufficiently small Aq as can be seen in 
Fig. El 

The general solution of Eq. is a difficult task, but 
simple analytic expressions can be found in the special 
cases of negligible and dominating atom-atom interac- 
tion. Omitting the last term in Eq. |^ it is straight- 
forward to see that |mcff| controls the magnitude of the 
dispersion term and thus the timescale of the wave packet 
broadening. A change in sign of rrics corresponds to time 



reversal of the dynamics in a frame moving with veloc- 
ity Vg. In the regime where the atom- atom interaction 
is dominating, e.g. during the initial expansion of a con- 
densate, the evolution of the envelope function can be 
found in standard nonlinear optics textbooks [l^ and in 
form of scaling solutions in the context of Bose-Einstein 
condensates [I3 • Note that in this regime the kinetic en- 
ergy term is still relevant and thus a change of the sign 
of the effective mass will reverse the dynamics. 

In the following we will analyze the obtained experi- 
mental results in more detail. The initial wave packet 
shown in Fig.^i) is characterized by Axq = 14.8(6) /im 
{Ax is the r.m.s width of a Gaussian fit) and shows 
the density distribution of a Bose-Einstein condensate 
of 2 ■ 10* atoms in a three-dimensional dipole trap. Be- 
fore releasing the atomic cloud into the one dimensional 
waveguide, a weak periodic potential along the waveg- 
uide is adiabatically ramped up to Vq — 2.8(2) Eq within 
6 ms. This turns the initial Gaussian momentum dis- 
tribution of the atoms into a Gaussian distribution of 
quasimomenta w{q) centered at gc = in the lowest en- 
ergy band of the lattice with a corresponding effective 
mass TOoff = 1.25 TO. The density distribution shown in 
Fig. ^p) is a result of propagation within the stationary 
periodic potential for i = 26 ms and exhibits a spread 
of 5x :— ^ Aa;(i)2 — Axg = 18.4(12)/^m in contrast to 
8xf = 20.2(14) /im for expansion without periodic po- 
tential. The resulting ratio 5xf/5x = 1.10(15) indicates 
that the evolution is dominated by the nonlinearity, in 
which case one expects Sxf/6x « ^Jrn^^Jm, = 1.11 jl4l |. 
In the case of linear propagation one expects 5xf/5x = 
TOos/to = 1.25. 

The dynamics in the anomalous dispersion regime are 
investigated by initially accelerating the periodic poten- 
tial within 3 ms to a velocity v = vg '■— hG/2m, thus 
preparing the atomic wave packet at the edge of the 
Brillouin zone, where TOeg = —0.5 to. The velocity is 
kept constant during the subsequent expansion. In the 
regime of negative mass a condensate exhibits collapse 
dynamics similar to condensates with attractive atom- 
atom interaction. Two-dimensional calculations for our 
experimental situation reveal that this collapse happens 
within the initial 3-6 ms of propagation. This leads to a 
fast reduction of the density and therefore of the nonlin- 
earity due to excitation of transverse states. An indica- 
tion of the population of transverse states is the observed 
increase of the transverse spatial extension of the wave 
packets for anomalous dispersion by almost a factor of 
two. The optical resolution of our setup does not allow 
for a quantitative analysis of the transverse broadening. 
The measured spread of 5x = 38.5(15) /xm after 23 ms of 
expansion leads to a ratio 5xf/Sx = 0.46(5). This ra- 
tio represents an upper limit due to the initial collapse, 
indicating that a coherent nonlinear evolution in one di- 
mension cannot describe the experimental results. 

In the case of dispersion management Fig. QJi) the 
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FIG. 2: (a) Band structure in the first Brillouin zone 
for atoms in an optical lattice with Vb = 1.2 Eg (solid), 
parabolic approximation to the lowest energy band at g = 
G/2 (dashed), corresponding group velocity (b) and effective 
mass (c) in the lowest energy band. The vertical dashed lines 
at g = g^ indicate where |mcff| — oo. (d-f) Spatial densi- 
ties of the wave packet after t = 60 ms of propagation with 
Qc = G/2 for different Vb- The position x along the waveguide 
is measured in the moving frame of the optical lattice. The 
solid lines represent the theoretical predictions using linear 
propagation with the exact band structure and the quasimo- 
mentum distribution given in graph (c). The dashed lines in 
Graph (d-f) represent the prediction of the constant effective 
mass approximation. 



wave packet was first subjected to normal dispersion 
for 17 ms at = 0. The time of subsequent prop- 
agation with anomalous dispersion at qc = G/2 was 
adjusted to achieve the minimal wave packet size of 
Ax = 15.4(2) /im. The minimum was achieved for times 
ranging from 7 ms to 9 ms which is in rough agreement 
with the expected time resulting from effective mass con- 
siderations -^^0.5/1.25 X 17 ms — 10.7 ms. The experi- 
mental results exhibit a rapid increase of the wave packet 
width for overall propagation times longer than 27 ms 
also indicating collapse dynamics. 

Before presenting the detailed experimental studies of 
the applicability of the constant effective mass approxi- 
mation and evidence of infinite masses, we turn to a brief 
description of our setup. 

The wave packets have been realized with a *^Rb 
Bose-Einstein condensate. The atoms are collected in 
a magneto-optical trap and subsequently loaded into a 
magnetic time-orbiting potential trap. By evaporative 
cooling we produce a cold atomic cloud which is then 
transferred into an optical dipole trap realized by two 
focused Nd:YAG laser beams with 60 /xm waist crossing 
at the center of the magnetic trap. Further evaporative 



cooling is achieved by lowering the optical potential lead- 
ing to pure Bose-Einstein condensates with up to 3 • 10^ 
atoms in the \F = 2,mp = +2) state. By switching off 
one dipole trap beam the atomic matter wave is released 
into a trap acting as a one-dimensional waveguide with 
radial trapping frequency luj_ = 27r • 80 Hz and longitudi- 
nal trapping frequency o^n — 2tt ■ 1.5 Hz. 

The periodic potential is realized by a far off-resonant 
standing light wave with a single beam peak intensity 
of up to 5W/cm^. The chosen detuning of 2nm to the 
blue off the D2 line leads to a spontaneous emission rate 
below 1 Hz. The frequency and phase of the individual 
laser beams are controlled by acousto-optic modulators 
driven by a two channel arbitrary waveform generator 
allowing for full control of the velocity and amplitude 
of the periodic potential. The light intensity and thus 
the absolute value of the potential depth was calibrated 
independently by analyzing results on Bragg scattering 
01 and Landau Zener tunneling 01 . 

The wave packet evolution inside the combined poten- 
tial of the waveguide and the lattice is studied by taking 
absorption images of the atomic density distribution af- 
ter a variable time delay. The density profiles along the 
waveguide, n{x,t), are obtained by integrating the ab- 
sorption images over the transverse dimension z. 

Since the assumption of a constant effective mass used 
so far is only an approximation, it is important to check 
its applicability in experiments. Therefore we investi- 
gate the dynamics of wave packets prepared at the Bril- 
louin zone edge {qc ~ G/2) for different potential depths. 
The observed density profiles after 60 ms of propagation 
are shown in Fig. Eli-f). While both the initial wave 
packet shape n{x, 0) and the quasimomentum distribu- 
tion w{q) are measured to be approximately Gaussian, 
the wave packet changes its shape during evolution. We 
attribute this distortion to the invalidity of the constant 
effective mass approximation, which assumes that the 
populated quasimomenta experience the same negative 
curvature of Eo{q). Since the range of quasimomenta 
fulfilling this criterion becomes smaller with decreasing 
modulation depth, a more pronounced distortion of the 
wave packet shape for weak potentials is expected (see 

Fig.m). 

This explanation is confirmed more quantitatively by 
comparing the observed wave packets with the results of 
an integration of Eq. |^ neglecting the nonlinear term. 
Since the initial collapse of the condensate cannot be de- 
scribed by a linear theory, we take a Gaussian function 
fitted to the density distribution measured at 20 ms as 
the initial wave packet for the numerical propagation. 
Due to the fact that this is not a minimum uncertainty 
wave packet we add a quadratic phase in real space such 
that the Fourier transform of the wave packet is consis- 
tent with the measured momentum distribution. In first 
approximation this takes into account the initial expan- 
sion including the repulsive atom-atom interaction. For 
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FIG. 3: Group velocities of steep edges emerging from an 
initial wave packet with significant overlap with in quasi- 
momentum space. The measured velocities of the indicated 
positions (arrows in insets) agree very well with the expected 
maximum and minimum velocity in the lowest band (solid 
line) corresponding to the infinite masses. The dashed line 
represents the prediction of the weak potential approxima- 
tion 9]. For potentials smaller than 2 Eg (diamonds) data 
arc obtained by preparing the initial wave packet at g = G/2 
leading to two steep edges (see inset). For higher potentials 
(circles) the wave packet is prepared at q — G/4 to ensure 
population of the quasimomentum corresponding to infinite 
mass. 



the subsequent propagation of 40 ms we use the full ex- 
pression for E(){q) which is obtained numerically. Finally 
we convolute the obtained density distribution with the 
optical resolution of our setup. In Fig.|3i-f) we compare 
the data with the linear theory described above (solid 
line) and with the constant effective mass approximation 
(dashed line). Clearly the effective mass approximation 
cannot explain the observed distortion and it strongly 
overestimates the expansion velocity for weak potentials. 
Additionally, for small potential modulation depths new 
features appear in the central part of the wave packet 
which cannot be explained using the linear theory. We 
are currently investigating these features in more detail. 

The observed distortion is mainly a consequence of 
another very interesting feature of the band structure: 
the existence of |TOcff| — oo for certain quasimomenta 
q = (see Fig. Et) . A diverging mass implies that the 
group velocity is extremal and the dispersion vanishes as 
can be seen from Eq. As a consequence an atomic 
ensemble whose quasimomentum distribution is overlap- 
ping q — q^ will develop steep edges as can be seen 
in Fig. [51i) . These edges propagate with the maximum 
group velocity of the lowest band. 

The systematic investigation of the velocities of the 
edges is shown in Fig.Ofor different values of Vq. In order 
to get a significant overlap of w{q) with q = g^, we pre- 



pare atomic ensembles with Aq = 0.17 G/2 at qc = G/2 
realized by Bose-Einstein condensates of 2 • lO** atoms 
with a spatial extension of Ax = 15 /im. The veloci- 
ties of the edges are derived from two images taken at 
20 ms and 60 ms, respectively. In each image the posi- 
tion of the edge is evaluated at the levels indicated by 
the arrows in the insets of Fig. O (50% and 25% of the 
maximum density). Since the momentum spread is too 
small to populate the infinite mass points for potentials 
deeper than 2Eg the atomic ensemble was then prepared 
at qc = G/4 by accelerating the periodic potential to the 
corresponding velocity v = ug/2. The resulting wave 
packet shapes are asymmetric exhibiting a steep edge on 
one side which becomes less pronounced for potentials 
deeper than 5Eg. The obtained experimental results in 
Fig. 131 are in excellent agreement with the numerically 
calculated band structure predictions. In contrast to the 
good agreement of the maximum velocity for all poten- 
tial depths we find that for Vq > 5 Eq the group velocity 
of the center of mass is only 10% of the expected veloc- 
ity deduced from the dispersion relation at qc = G/4. 
This could be an indication of entering the tight bind- 
ing regime where the nonlinear effect of self trapping, 
i.e. stopping and non-spreading wave packets, has been 
predicted by Trombettoni et al. 0. We are currently 
investigating the transport properties in this regime in 
more detail. 

In conclusion, we have demonstrated experimentally 
that the dispersion of atomic matter waves in a wave- 
guide can be controlled using a weak periodic potential. 
Matter wave packets with positive, negative and infinite 
effective masses are studied in the regime of weak and in- 
termediate potential heights. The preparation of matter 
waves with engineered dispersion (rricff < 0) is an im- 
portant prerequisite for the experimental investigation 
of atomic gap solitons and other effects arising from the 
coherent interplay of nonlinearity and dispersion in peri- 
odic potentials. 
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